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We study phase transitions and the nature of order in a class of classical generalized 0{N) 
nonlinear cr-models (NLS) constructed by minimally coupling pure NLS with additional degrees of 
freedom in the form of (i) Ising ferromagnetic spins, (ii) an advective Stokesian velocity and (iii) 
multiplicative noises. In examples (i) and (ii), and also (iii) with the associated multiplicative noise 
being not sufficiently long-ranged, we show that the models may display a class of unusual phase 
transitions between stiff and soft phases, where the effective spin stiffness, respectively, diverges and 
vanishes in the long wavelength limit at two dimensions (2d), unlike in pure NLS. In the stiff phase, 
in the thermodynamic limit the variance of the transverse spin (or, the Goldstone mode) fluctuations 
are found to scale with the system size L in 2d as InlnL with a model-dependent amplitude, that 
is markedly weaker than the well-known In L-dependence of the variance of the broken symmetry 
modes in models that display quasi-long range order in 2d. Equivalently, for A = 2 at 2d the 
equal-time spin-spin correlations decay in powers of inverse logarithm of the spatial separation with 
model-dependent exponents. These transitions are controlled by the model parameters those couple 
the 0{N) spins with the additional variables. In the presence of long-range noises in example (iii), 
true long-range order may set in 2d, depending upon the specihc details of the underlying dynamics. 

Our results should be useful in understanding phase transitions in equilibrium and nonequilibrium 
low-dimensional systems with continuous symmetries in general. 


I. INTRODUCTION 


Understanding phase transitions and order in low di¬ 
mensional systems has remained a topic of considerable 
theoretical interests. For two-dimensional {2d) systems in 
equilibrium with continuous symmetries and short-range 
interactions, there is no phase transition to a low temper¬ 
ature (T) broken symmetry phase with long range order 
(LRO) at any T 0. This is a consequence of the well- 
known Mermin-Wagner theorem (MWT) [I|. As a result 
in the thermodynamic limit (TL), there is no statistically 
flat phase of a fluid membrane or a 2d crystal in thermal 
equilibrium. In contrast, long-range order in 2d equilib¬ 
rium systems can exist at a finite T, if there are effective 
long-range interactions present in the system which can 
take the system out of the validity of the MWT. No¬ 
table examples include tethered membranes 0, where 
one finds a finite-T crumpling transition and a low-T 
statistically flat phase with long-range orientational or¬ 
der 0 at 2d. Symmetric inhomogeneous fluid membranes 
form a closely related example, see Ref. Q- Here a phase 
transition is obtained between a phase with a finite per¬ 
sistence length and a phase with a diverging persistence 
length (or equivalently, a stiffening transition), tuned by 
the coupling constant that couples the local curvature 
with the local inhomogeneity order parameter. However, 
true LRO does not exist even in the phase with diverging 
persistence length in Ref. The latter phase is found 
to be characterized by a diverging bending modulus in 
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the long wavelength limit at 2d. 

The classical 0{N) nonlinear cr-model Q (hereafter 
NLS) has found wide usage in a variety of topics, in¬ 
cluding condensed matter systems @ , mathematical 
physics @ , particle physics Q and cosmology Q . This is 
a paradigmatic model with continuous symmetries 0,0, 
in which the consequence of the MWT is succinctly vis¬ 
ible in a low-T expansion. At d = 2 in equilibrium, one 
finds 0 that the model has no finite-T ordered state with 
LRO. In this article, we ask: what can introduce phase 
transitions and order in NLS at 2dl Taking cue from 
the examples of tethered membranes, we expect that the 
predictions of MWT may be bypassed in the presence of 
long-range interactions between the 0{N) spins. Since 
pure NLS has only short-range interactions between the 
spins, additional degrees of freedom which may create ef¬ 
fective long-range interactions, are required to invalidate 
the prediction of the MWT here. These additional fields 
can be interpreted as representing spatio-temporally non¬ 
trivial environments coupled to NLS. The model parame¬ 
ters defining pure NLS now get modified or renormalized 
by the coupled additional degrees of freedom. Clearly, 
any putative LRO in a 2d generalized NLS should depend 
upon the specific nature of the additional degrees of free¬ 
dom and their couplings to the 0{N) spins, i.e., the re¬ 
sults should be mo del-dependent. Since we are interested 
in a question of general principle, it is useful to study sim¬ 
ple reduced models, where explicit calculations may be 
done in straightforward ways. To this end, we generalize 
NLS and construct three different simple variants of NLS 
by coupling it to additional fields through generic local in¬ 
teractions, one in equilibrium and the others out of equi¬ 
librium: We consider a system of NLS, that is (i) thermo¬ 
dynamically coupled with Ising spins - hereafter Model 
I, (ii) dynamically coupled with a Stokesian velocity field 
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that is affected by the feedback from the 0{N) spins - 
hereafter Model II, and (hi) dynamically coupled with 
multiplicative noises through symmetry-allowed minimal 
coupling - hereafter Model III. For Model I and Model II, 
we consider only nonconserved dynamics while for Model 
III we consider both nonconserved and conserved dynam¬ 
ics of NLS. In a low noise variance expansion that gen¬ 
eralizes the low-T expansion for the equilibrium NLS, we 
illustrate the possibility of order in the NLS with ad¬ 
ditional fields at 2d. Our principal results are (i) both 
Model I and Model II and nonconserved Model III with¬ 
out any long range noise have the lower critical dimension 
dL = 2; for d > di all these have a critical point separat¬ 
ing a “high noise” paramagmetic phase with short-range 
order (SRO) and a “low noise” ferromagnetic phase with 
LRO, (ii) in 2d, both Model I and Model II, and the non¬ 
conserved version of Model III (without any long-range 
noise) with appropriate choices of the model parameters, 
show phase transitions, tuned by the coupling constants, 
from a disordered soft phase with SRO, where the ef¬ 
fective spin stiffness of the 0{N) spins vanishes over a 
finite system size to a stiff phase, where it stiffens sig¬ 
nificantly, diverging as InL for a system of linear size L. 
As a result, in TL the variance of the fluctuations of the 
transverse components of the spins depend on L as In In L 
at 2d, a dependence weaker than the standard InL de¬ 
pendence of the variance of the elastic degree of freedom 
in an elastic Hamiltonian that displays quasi long-range 
order (QLRO) in 2d. For = 2 at 2d this implies a spa¬ 
tial decay of the equal-time spin correlations in powers 
of the logarithm of the spatial separation with model- 
dependent exponents. The model-dependence of these 
exponents are reminiscent of the model-dependent expo¬ 
nents that characterizes the well-known algebraic decay 
in QLRO; nevertheless, the spatial decay here is markedly 
slower than the algebraic decay in QLRO. In contrast, the 
conserved version of Model III without any long-range 
multiplicative noise admits no stiff phase at 2d. (ii) Non¬ 
conserved Model III with long-range multiplicative noises 
can display true long-range order at 2d, controlled by 
the noise amplitude and the relevant coupling constant. 
However, the conserved version of Model HI with long- 
range multiplicative noises does not display any LRO in 
2d. In addition, we also calculate the dynamic exponent 
z that characterizes the correlators of the transverse spin 
fluctuations in the stiff phases at 2d and at the unstable 
fixed point (FP) separating the paramagnetic and ferro¬ 
magnetic phases for d > 2. The remainder of this article 
is organized as follows. In Sec. m we provide a short 
review of NLS. In Sec. Ell 113 and 13 we construct and 
analyze our Model I, II, and HI, respectively. In Sec. IVII 
we summarize our results. Some of the technical details 
are available in Appendices at the end. 


II. SHORT REVIEW ON NLS 

The statistical mechanics of NLS with an N- 
component spin <I> = {(fi, in the presence of an 

external magnetic field hi, i = 1,N in d-dimensions is 
described by the free energy functional Ta 0 

= d''a;[«(V„$)2 - h,cj)i], (1) 

where k is the spin stiffness; we impose = 1. By means 
of a perturbative low-T expansion, it has been shown 0 
that the transition temperature Tq separating the low-T 
ferromagnetic phase and the high-T paramagnetic phase 
is given by 


at d = 2-l-e, e > 0. Thus, Tq vanishes at 2d (i.e., d^ = 2), 
precluding any finite T ordered phase with LRO. For 
N = 2 at 2d, To is indeterminate, suggesting that de¬ 
grees of freedom not included in NLS, e.g., topological 
defects and amplitude fluctuations, should destroy the 
order. Similar analyses have been done on the quantum 
version of the 0{N) nonlinear cr-model; see Ref. [^. We 
now ask whether or how these results may be significantly 
modified in the presence of additional degrees of freedom. 
We construct Models I, II and HI as above to address 
this issue. To set up the notations clearly, we denote the 
spin components of an 0{N) spin by Roman indices and 
the components of a vector or an operator in the real or 
Fourier space by Greek indices in the remaining part of 
this article. 


III. MODEL I: NLS COUPLED WITH ISING 
SPINS 

We study the equilibrium nonconserved relaxational 
dynamics of NLS coupled with Ising spins. In the spirit 
of the Landau-Ginzburg coarse-grained approach, we rep¬ 
resent the Ising spins by a scalar field ^/'(x), where x is 
the coordinate in a d-dimensional space. We start with 
the combined free energy functional T/ for NLS coupled 
with the Ising spins: Assuming a minimal symmetry- 
permitted coupling between the Ising spins and the 0{N) 
spins, we write 

Tz = y d''cr[|(V„$)" - h • $ + + A(V„$)V" 

+ (3) 

where, as before, <I> = ((/>i,..., ^at), is an A-component 
vector of unit modulus,i.e., = 1. Here, r = T — T^, T 

being the temperature, T^ is the mean field critical tem¬ 
perature for the second order Ising magnetic transition, 
coupling constant u > 0. Setting ijj = 0 yields Eq. (HD 
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above; k is the spin stiffness. For N = 2, © is related 
to a Ginzbug-Landau model recently proposed in stud¬ 
ies on high temperature superconductors (ll| . In ([3|) the 
coupling between (j)i and Tp are controlled by the coupling 
constant A; With T > Tc, for all A > 0, is minimized 
by uniform states of the fields, i.e., pi = const, and '0 = 0, 
where as, with large enough A < 0, such a uniform state 
gets thermodynamically unstable. We choose A > 0 here. 
It is constructed such that Eq. © is invariant under rota¬ 
tion of the 0{N) spin pi in the spin space and inversion of 
the Ising spins p —>■ —p. Usual order-disorder transitions 
and relaxational dynamical critical behavior near critical 
points in a related model in higher dimensions (e.g., three 
dimensions) have been studied in Ref. [H , where a va¬ 
riety of temperature driven phase transitions separating 
the disordered and the ordered phases are discussed. Un¬ 
like Ref. [l^, we here assume that the 0{N) spins are 
in the ordered state (and hence put a fixed length con¬ 
straint on the 0{N) spins) and study the effects of the 
Ising spins coupled with the 0{N) spins via generic sym¬ 
metry allowed couplings on the putative ordered state of 
the 0{N) spins at 2d. Since the Ising spins undergo a 
usual paramagnetic to ferromagnetic transition at T^, we 
focus on the state of the assumed order of the 0{N) spins 
B.t T = Tq. 

We are interested in studying the purely relaxational 
nonconserved dynamics. The equations of motion for $ 
and p are respectively given by: 

^ = -r^+0 = r [kvH + 2AV;3(0"V;3$) + h] +e, 

(4) 

and 




-r. 


STi 

-rp - 2XP{V+ VV + 


T]iP>) 


Here, F and r 2 are the kinetic coefficients of $ and p, 
respectively. A quick inspection of Eq (jT]) reveals that p- 
fluctuations contribute positively to n through the term 
with coefficient A. Whether or not this can compete with 
the well-known thermal softening of n that is present even 
in pure NLS Q can only be determined after a thorough 
calculation. This remains a major goal of this work. 
Stochastic functions 9 and r] are zero mean Gaussian 
white noises with variances given by [we set Boltzmann 
constant fcs = 1 ] 


(0,(x,t)0,(O,O)) = 2DT6,j5{^)5{t), (6) 

(7?(x,t)7?(0,0)) = 2T2T5{yi)5{t). (7) 


We are considering equilibrium dynamics of the cou¬ 
pled system. Hence, the system obeys the Fluctuation- 
Dissipation Theorem (FDT) BIS' This enforces D — 
T. Noting that T/k is a dimensionless number at 2d, we 
intend to calculate the effective or renormalized k in a 
low temperature expansion {D/k <C 1) up to one-loop 
order. For simplicity we treat the fluctuations of p up to 


the harmonic order, i.e., we set u = 0 [l3|. For simplicity, 
we now set r 2 = 1. 

It is clear that a uniform state of the 0{N) and Ising 
spins, where the 0{N) and Ising spins are perfectly 
aligned among themselves, minimizes J-j. Such a state 
may be parameterized as $ = (1,0, ...,0). Thermal fluc¬ 
tuations should reduce the component of $ along the 
direction of order. Now assume an ordered state at low- 
T where the 0{N) spin system is assumed to be ordered 
in a given direction. A convenient parameterization of 
$ for such a state is $ = (cr, tt); thus, with cr « 1 
and TTi,i = 1,..,A^ — 1 are the small transverse fluctu¬ 
ations; hence, tt is an — 1 component vector. We write 
a = \/l — TT^ Ri 1 — where tt^. We impose 

the condition = 1 on the dynamics via imposing a 
(5-function, 5i = i5 [cr^ -I- tt^ — l] on the generating func¬ 
tional to be constructed (see iDelow) from (j?]) and ([S]). 
Notice that the fixed length constraint = 1 essentially 
implies $ • dt^ = 0, i.e., the variation in $ takes place 
normal to $, or, in other words, there are (A^ — 1) inde¬ 
pendent transverse fluctuating degrees of freedom. The 
latter are to be associated with defined above with 
fj « 1. Thus, the stochastic force 9 = dt^ + T5iFi/5^ 
in o must also be purely transverse. This is ensured by 
imposing a second i5-function (52 = i5[$ • d>] on the associ¬ 
ated generating functional MM, where $ = {d, tt) is 
the dynamic conjugate of [l5|. Clearly, contributes 
only when a = \/l — tt^. Further, with $• d) = dd + Tf -tv, 

here di is the 


\/l —71 


82 contributes for d = — 
dynamic conjugate of tt^. 

After inclusion of the constraints in the form of the 
(5-functions 81 and 82 , the generating functional for the 
system is given by (assume h is aligned along a with a 
magnitude hi) 


Zj = J 'Da'Dd'D'K'Dd'Dp'Dp exp [S'/] 

= / 'Da'Dd'Dn'Dtr'Dp'Dp exp[ / d^xdtdd 


D 


D 


+ J d‘^xdt—TtjTt_ 

/ r 1 d(T 

d^xdtd < “ kV^ct — 2W/jipPVpa) — hi 

— J d'^xdtdj — KV^TTj — 2XWp{p^'X/pTTj)'^ 

+ J d‘^xdtTpp — J d^xdtp{dtp + rp— V‘^p 

+ 2Xp{Vpa)^+2XP{Vp-K)^}] 


X 8 [a'^ + — l] 8 [da + diTTi\ 


( 8 ) 


The functional integrals over a and d in Zj above may 
be evaluated by using the two (5-functions, which yield 
two Jacobians. Overall the total Jacobian J evaluates to 
(1 — TT^), appearing at every x and t. Thus, the contri¬ 
bution of the total Jacobian to Z is a product of all these 
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terms,wz. fl 5 l|. 

- TT^) = exp[- ^ /dt lll(l - TT^)] 

X 

= exp[—p J df^xdt ln(l — (9) 

where, p is a number equal to the number of degrees 
of freedom per unit volume that introduced while going 
from summation to integration. We now expand the non¬ 
linear terms up to 0{T/k), since we are interested in a 
low-T one-loop expansion. We obtain the action func¬ 
tional after truncating to the linear order in (T/ k) 

rt f id iiT^ ^ ^ .. r ^ d^Tj t—,2 

Si = J d xdti — TTjTTj - TTj{-— - kV TTj 

- 2X\//sillily pTTj)} + ^{TTjTTjf 
1 1/^2 

+ - 2hi{l + 7r^/2)} -f 

— -f -I-2A^/>(Vp7r)^}-I-pTT^]; (10) 

see Appendix El for details. Notice that while the contri¬ 
bution from J nominally yields a quadratic term of the 
form pTT^, it is 0(T/k) [T^; see also Appendix lAl It is 
conyenient to work in the Fourier space; we define q and 
ui as the Fourier wavevector and frequency, respectively. 
At T > Tc, the fluctuations in ^ are short lived and 
small; as a result, the contributions to the measurable 
physical quantities from the ^-fluctuations are small for 
T > Tf, [ 13 . In contrast, the ^-fluctuations diverge in the 
long wavelength limit as T T^. Similarly, being a 
broken symmetry mode (i.e., a Goldstone mode), correla¬ 
tion (|7ri(q, t)p) diverges in the long wavelength limit at 
all T. Thus, contributions to the measurable quantities 
from '(/j-fluctuations can compete with those originating 
from the broken symmetry mode fluctuations only near 
the critical point, i.e., as T —>■ T^. Hence, in our subse¬ 
quent analysis below in this Section, we set D = T = Tc- 
In order to deal with these long wavelength divergences in 
a systematic manner, we employ Wilson momentum shell 
dynamic renormalization group (DRG) [Hi ; see also 
Ref. for detailed discussions on DRG applications to 
dynamic critical phenomena. To this end, we first in¬ 
tegrate out fields 7ri(q,w),^(q,w),7ri(q, w),'i/'(q, w) with 
wavevector A/6 < g < A, 6 > 1, perturbatively up to the 
one-loop order in (1101) . Here, A is an upper cut off for 
wavevector. This allows us to obtain the ’’new” model 
parameters corresponding to a modified action Sf with 
an upper cutoff A/b < A. We obtain 


K< 

= K[l + A-k2-A], 

(11) 

r 

= ^[l + A], 

(12) 

1 

f< 

= ^[1 + A], 

(13) 

h< 

N -1 

= hi[l + ^A]. 

(14) 


Here, superscript ^ refers to the parameters in the action 
with a reduced upper cut off A/b. Notice that (D/F)< 
and l/F^ have the same form, which is a consequence of 
the FDT. Furthermore, 

d<^q D ~ d‘^q T 

J A/b {i^q^ + hiY J (2ttY r + 

(15) 

where we have used the forms of the correlators given in 
(IbI) . Glearly, at 2d, A is log divergent for small hi while 
A is log divergent at r(= T — Tc) = 0. 

In order to extract the renormalized parameters, we 
then rescale wavevectors and frequencies according to 
q' = 6q and lo' = b^uj, where 2 is the dynamic exponent. 
Under these rescalings, fields (/iiYi also scale. We write 
4>i —t iYi under the above rescalings. FDT 

in the present model implies that Im{TTiTTi) and uj{TTiTTi) 
must scale in the same way under the above rescalings 
{Im refers to the imaginary part). This consideration 
yields ^ = ^b~^. With this, rescaling factors for the var¬ 
ious model parameters may be obtained. These lead to 
the renormalized parameters 


h'l 



U-d-: 


ahi 


1+(^A 


<^6 


— (d+2+2) 


. A 2A;r' 

1 + A + —A 

K 


[1 -(-A]. 


(16) 

(17) 

(18) 
(19) 


Notice that with i = ^b the left hand sides of Eqs. dUD 
and (fT^ scale in the same way. This is a consequence of 
the FDT as mentioned above. 

In general, the external magnetic field h itself is an 
0{N) vector that couples with the 0{N) spin $ through 
the coupling with the external magnetic field given by 
— f hj(j)jd'^x in the free energy. Now consider the cou¬ 
pling — J hj'Kjd'^x in the free energy ([3]) which would 
produce a term f hjTTjd‘^xdt in the corresponding ac¬ 
tion functional m- Requiring this ’’external part” 
of the action to be invariant under rescaling implies 
J hj-iTj d‘^xdt = f h/n'j d’^x'dt'. Since the 0[N) symme¬ 
try of the problem with hi = 0 ensures that tt^ scales the 
same way as cr, we have 6,' = |/ii; see Ref. 0 for similar 
arguments for pure NTS.. Due to the 0{N) symmetry of 
the model for hi = 0, all components of hi must scale in 
the same manner. Hence, 


h'i = ihi. 


( 20 ) 


Using Eq. (HSl) and Eq. (uni, we now evaluate ^ and 


C = 6^^+^ 
i = b‘^ 


1 - 

.-A, 


( 21 ) 

( 22 ) 


D 

Y 
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These yield 


k' = nb'^ ^ 


2 X 

1 - (A^- 2)Ah-A 

K 


(23) 


We set D = T = Tc and are interested in the fluctuation 
corrections at d = 2 + e, e > 0. At the leading order in 
e, we evaluate A and A above at d = 2. We thus find 
A = ^^ln6 for small hi and A = |jln6. 

We now obtain the continuum flow equations for k and 
r. Let b = e\ 1^0 and d = 2+e. Using these definitions 
and Eq. (|23)l . the differential flow equation for k becomes 


dK \ (N- 2)Tc 2XTc 

— = K e -1- 

dl 2 ttk 2ttk 


(24) 


First consider N > 2. For N — 2 > 2A, there is a DRG 
fixed point (FP) given by dn/dl — 0. This yields at the 
unstable FP 


Tc 27re 

~K ^ (A - 2) - 2A’ 


(25) 


at d = 2 + e. Compare (l2^ with the corresponding result 
for pure NLS for the (reduced) transition temperature 
given by Thus, (1251) implies a (reduced) transition 
temperature T = Tc/ k. having a value T* ^ 0(e) at the 
FP that separates a low temperature ferromagnetic phase 
(ordered state), where the majority of spins are aligned 
and a high temperature paramagnetic phase without any 
alignment of spins at d = 2 + e. Since T* = 0 at 2d 
(e = 0), we have d^ = 2. This physical picture holds 
for A > 2 + 2A for which T* drops to zero at 2d. At 
A = 2 + 2A, T* is indeterminate at 2d. This is the 
analog of the behavior of the transition temperature for 
pure NLS with A = 2. It is understood that topological 
defects, e.g., vortices and amplitude fluctuations destroy 
the order for A = 2, d = 2 in pure NLS Since A must 
be an integer, A may be increased in steps of unity, where 
as A is any real number. Thus, A = 1/2 yields A = 3 
such that at d = 2, T* is indeterminate for Model I with 
A = 3. It will be interesting to study what role, if any, 
topological defects play to destroy the order there. 

In the above analysis, we have assumed the existence 
of an unstable FP, which may be ensured by choosing A, 
a free parameter in the model, appropriately. For A < 
2 + 2A, dn/dl > 0 thus precluding any FP. To proceed 
further, we define a critical Ac(e), at which ^ vanishes, 
as a function of the e parameter given by 


Ac(e) 


(A — 2) TTKC 
2 


(26) 


Defining AA(e) = A —Ac(e), (IMl) at d = 2 may be written 
as 


dK TcAX{e = 0 ) 

dl TT 


order phase transition; nonetheless, there are significant 
differences that we elaborate below. 

Consider now the scale-dependent or renormalized k{1) 
at an off-critical point (A ^ Ac or A A ^ 0). Clearly from 
Eq. (l27l) . for AA(e = 0) < 0, k{1) reduces linearly with I 
at 2d and ultimately becomes zero at a particular length 
scale depending on kq and A. Thus even for NLS coupled 
to Ising spins at Tc, with A < Ac(e = 0), k vanishes for 
a sufficiently large system at 2d. The discrete recursion 
relation for k at 2d is given by 

n = Kii- (^~^)^^ in6 + ^ln6, (28) 

ZTT TT 

where kq is a microscopic stiffness. Further writing ln5 = 
ln(A/g') = —In(gao), where uq a microscopic length, we 
obtain the ^-dependence of k : 

«(?) = Ko + ^ln(gao) - ^^ln(gao). (29) 

ZTT TT 

Evidently, the second and third terms of (|2^ reveal the 
competition between thermal softening of k and its stiff¬ 
ening due to the '(/-fluctuations. Not surprisingly, the 
wavevector-dependent effective or renormalized k does 
not depend upon P, a kinetic coefficient. This is again 
a consequence of the FDT-obeying dynamics. Equa¬ 
tion (pn)) allows us to define a persistence length ( ^ 1/q, 
such that k(() = 0 (see, e.g., [^[l3 for similar persistence 
lengths in different problems). This yields 


C = aoexp 


27rKo 

(A - 2)Tc - 2XTc 


(30) 


Hence, as A —>■ Ac_ (e = 0) = in 2d, the persistence 

length, C —>■ oo. Civen the physical interpretation that 
C is a typical length scale at which k{C) = 0, the system 
remains ordered over a length scale smaller than (/. Thus 
for system sizes larger than (/, there will only be SRO. 
We also define persistence length ((NTS) for pure NLS 
(A = 0) 


((NTS) = aoexp 


27rKo 

(A - 2)Tc 


(31) 


Clearly, even for A < Ac, ((NTS) < (. Thus, even in the 
soft phase, the persistence length is larger than that for 
pure NLS, implying that A in general favors order. 

In contrast for A > Ac(e = 0), k grows on successive 
applications of momentum shell DRC. Thus, at larger 
scales NLS is supposed to appear more ordered than what 
it is at smaller scales. Clearly at higher dimensions, as 
e increases, Ac(e) decreases. Hence at higher dimensions, 
weaker 0(A)-Ising couplings are enough for k(1) to grow. 
Consider now the explicit form for K(q) in the limit q —>■ 0 
for A > Ac(e = 0) at 2d. For A — Ac(e = 0) = AA(e = 
0) > 0 at 2d, we have. 


giving an unstable FP A = Ac(e = 0) = (A — 2)/2 at 
2d. This is reminiscent of a usual equilibrium second 


/ \ dcAA(e = 0), , , 

K(q) =-^ln(gao) + Kq, (32) 

TT 





















6 


which diverges logarithmically in TL (g —^ 0). In order 
to comment on the nature of order, if any, displayed by 
Model I, we look at the variance of the transverse spin 
fluctuations tt, given by Aq = (|7ri(q,t)p), now defined 
in terms of the renormalized bending modulus K{q). We 

f (Pq Tc 1 f dq 

J {2ttP K{q)q^ + hi 2AA(e = 0) J q\n{qao) 

m(bo)i“ii“(“"/-^)i = 2 aa(! = o) ‘"‘"^’ 

(33) 

in the limit of large systems, i.e., L ^ oq, where L is a 
lower momentum cut-off ^ linear system size. Regardless 
of the value of Tc, (155)) holds as long as A > Ac(e = 
0). Thus similar to Ref. Q, Aq diverges as L —>■ oo, 
but very slowly and remains finite for any finite value 
of L. In contrast for pure NLS, K{q) vanishes for a low 
enough q. Hence for the pure system, Aq will diverge 
even for a finite system size. For d > 2, it may be shown 
straightforwardly that Aq is finite in TL. Thus, true LRO 
exists for d > 2. A phase diagram of the system in the 
A—e plane is shown in Fig.[l] with the line of FPs given by 
A = Ac(e), SRO for A < Ac(e), LRO for A > Ac(e), e > 0 
and stiff phase for A > Ac(e = 0) . 


find 
Aq = 



FIG. 1: Phase diagram in the e —A plane for N = 2 (Ac = 0). 
The stiff phase exists along the A-axis (broken line, e = 0, A > 
0). Symbol O which marks the origin (0,0), corresponds to 
pure NLS with SRO at 2d. 


It is insightful to consider the flow of the reduced tem¬ 
perature T = TcIk av 2d. We find 


df f2.N-2 2A 
dl 2 tt 2tt 


(34) 


For A = 0 and iV > 2, T = 0 is the only DRG FP. 
This is, however, an unstable FP (expected), such that 
even for an arbitrarily small (microscopic) T = Tq, scale- 
dependent T{1) grows. Integrating (1541) 

— Y ~ 


where Aq = and C is a constant of integration. 

Assume Aq > 0. Set at Z = 0, l/T = 1/Tq > 0, we have 


_l=AoZ-i. 

T To 


(36) 


Thus, as I —>■ I/(AoTo) In^ from below, T —>■ oo. This 
implies that for any microscopic Tq > 0, the large scale 
properties of the system is identical to that of a system at 
T —>■ oo, indicating disordered phase at all T > 0. Since 
Tc for an Ising system at 2d is always larger than zero 
(there is a finite temperature Ising ferromagnetic phase), 
our Model I at T = T^ is always in its disordered phase 
for A = 0. This scenario changes for Aq < 0, or, A > Ac 
at 2d. Flow equation (1551) . upon integration, now yields 

i = IA|l + A (37) 

Evidently, for Z —?> oo, T —^ 0. Therefore, for A > Ac(e = 
0), scale-dependent T(Z) flows to zero as Z —>• oo, where 
as for A < Ac(e = 0), T(Z) flows to infinity as Z —>■ In^. 
Thus for A > Ac(e = 0), the properties of the system at 
the largest (formally TL) scale is expected to be same as 
that of the system at zero temperature. This suggests 
the existence of order in the system. On the other hand, 
for A < Ac(e = 0), the behavior of the system at large 
scale is same as that at infinite temperature, indicating 
a paramagnetic phase with SRO only. Thus, as A rises 
from 0 to a high value through Ac(e = 0), the system 
undergoes a phase transition from a disordered (para¬ 
magnetic) phase to an ordered phase through a critical 
point at A = Ac(e = 0). This is reminiscent of the usual 
temperature driven order-disorder second order transi¬ 
tion; however, in the presence case, T is kept fixed at Tc 
and A is the tuning parameter. 

We now discuss the nature of the order in the ordered 
phase, and whether it is identical to the magnetic (i.e., 
ferromagnetic) order in usual magnetic systems. The 
usual magnetic transition is conveniently described by 
the magnetic order parameter, which is non-zero in the 
ferromagnetic phase but zero in the paramagnetic phase. 
Since we assume order in the “I” direction in the spin 
space, we take m = (a) = {\/l — tt^) Ri I — (7r^)/2 as the 
order parameter. Clearly, in the disordered phase m Ri 0 
for a system with a linear size approaching Interest¬ 
ingly, even in the ordered phase, m r: 0 in TL, due to 
the divergences found in Aq above; see Eq. (1551) . Thus, 
m remains zero on both sides of the critical point given 
by A = Ac, and hence cannot work as an order parameter 
here [1^. Instead, we take 

O = [ln(e/ao)]-i = -TcAA/(^ko), (38) 

with A < Ac and 0 = 0 with A > Ac as the order param¬ 
eter at 2d. This is similar to the order parameter defined 
in the context of a crunmled to stiff phase of a heteroge¬ 
neous fluid membrane Q . Evidently, O in the soft phase 
rises smoothly from zero, as A is reduced from Ac. Thus, 
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with A as the control parameter, the order parameter ex¬ 
ponent is unity. In the stiff phase, O is naturally zero. A 
schematic plot of O versus A is shown in Fig. [2J Fig. [3] 
and Fig. 0] provide schematic phase diagrams of Model I 
in the N — X (with T = Tc) and X — T (for a fixed N > 2) 
planes at 2d [?7l| . 



FIG. 2: (color online)Schematic variation of the order param¬ 
eter O versus A in 2d for N > 2. 



N 


FIG. 3: (color online) Schematic phase diagram in A — A 
plane at 2d at T = Tc- The thick line (blue) corresponds to 
the equation Xc{e = 0) = {N — 2)/2. The stiff and soft phases 
are marked (see text). 


What happens for iV = 2 at 2dl In pure NFS, the 
transition temperature becomes indeterminate, suggest¬ 
ing the importance of topological defects and amplitude 
fluctuations in destroying LRO in the system. When the 
0 ( 2 ) spin is coupled to the Ising spin, the flow equation 


soft 


transition point 

1 

Stiff 

T • 

1 

' c • 


soft 


* 



_11S_*_ 

O A 


FIG. 4: (color online) Schematic phase diagram is A —T plane 
for A > 2 at 2d. Symbol, O = (0, 0) is the origin. Stiff phases 
exist along the horizontal thick line (red). Thick circle (blue) 
refers to the critical point (Ac,Tc). 


for k{1) reduces to 


dK XTc 

dl TT 


(39) 


giving K,{q) = kq — XTcln{aoq )/tt for all A > 0. Thus 
for any positive A, renormalized k{1) diverges in TL. As 
a result, Aq ^ [In lnL]/A at T = Tc and the 0(2) spins 
should appear more ordered than the XY model with 
QLRO at low T. In an equivalent picture, (|M)) for A = 2 
reduces to 


df 

m 



(40) 


Thus, T{1) always flows to zero for all A > 0. Hence, 
Model I with A = 2 is in its ordered stiff phase for all 
X. We now calculate the Debye-Waller factor and the 
equal-time spin-spin correlation Os(r) = (4>(x) • <l>(x')), 
for A = 2 in the stiff phase {T = Tc,X > Ac), where 
r = |x — x'l is assumed to be large. For A = 2, the 
transverse components tt has only one component, say, If; 
$ may be expressed as a complex number exp(fT') where 
ff'(x, t) is a real scalar field; we write = cosfl', $2 = 
sin'I' as the two components of $. For small transverse 
fluctuations in a nearly ordered state, associate 'F with tt 
(equivalently, parametrize $ = exp(j7f); see Appendix [Pl 
for formulation of the dynamics in terms of 4'). In the 
renormalized version of Model I, the variance of 7f(x, t) in 
the stiff phase (for all A > 0) in 2d is given by (1331) with 
e = 0 and A = 2. Now assuming “1” to be the ordering 
direction, the average order parameter 

(cos4') = exp[-i(ff(x,t)2)] = exp[-iAo] 

= exp[-^ln|lnL|] = exp(-IF), (41) 







in (renormalized) Model I for a system of (large) lin¬ 
ear size L. The Debye-Waller factor that determines the 
depreciation of the magnetization due to the noise in¬ 
duced fluctuations (thermal fluctuations for equilibrium 
systems) is given by exp(—2W), where 2W = Aq = 
(lnlnL)/(2A) in the limit of large L Q. Clearly, IT —>■ oo 
in TL here and hence the order parameter vanishes in TL. 
Note, however, that as A rises, W decreases for a fixed L, 
and thus, an increase in A reduces thermal depreciation 
of the magnetisation. Furthermore, 


Cs{r) = (cos[^'(x,t) - ^'(0,t)]) 

= Re{exp{i['S{x,t)-'if{0,t)])) =exp[-g(x)], 

where g{x.) = {[n{x,t) - n{0,t)]'^)/2 = 

f (|;;f7[(l^q(i)P)[l - exp(jq • x)], where we again 
have made the association between the phase field 
<i> and TT for small transverse fluctuations. Defining 
q • X = qr cos ip = u cos (^, we obtain in renormalized 
Model I 


5(x) 


I 

prA 

du ^ 

2A 

J rL~^ 

u \n{uao/r) 

I 

f 

du 

2A 

J rL ^ 

iiln(uao/r) 

I / 

^rA 

du 


2XJi itln(uao/r)’ 


J_ duJoju) 

2XJi u\n{uao/r) 

(42) 


where A ^ I/oq is an upper momentum cutoff; Jo{u) is 
the zeroth order Bessel function of first kind. Clearly, 
the first two terms in the last line of (1421) remain finite in 
TL, L —>■ oo. Thus in TL for r —>■ oo and neglecting the 
finite parts 

5 (x) = ^ln|ln(ao/r)|. (43) 

This gives 


1 

(Inr)^/^''' ’ 


(44) 


The last term vanishes at A < Ac(e) yielding 


dl 


kT = —kT 


2 -(- e — z — 


{N - 2)T, 


2ttk 


( 46 ) 


At the FP, this yields 2 = 2 -|- e — ■ When A < 

Ac(e), TJk = 2-Ke/\N - 2 - 2A] at the DRG FP, which 
gives z = 2 — (jv_^ 2 )- 2 A at the DRG FP. Thus z clearly 
depends on A and decreases as A rises for a given N and 
e. On the other hand, for A > Ac(e) and using (1^ . 
we find z = 2, since —?► 0 for L —>■ oo in the stiff 

phase. Notice that the dynamic exponent z^ for tp with 
a nonconserved relaxational dynamics is z^ = 2 to the 
linear order in u. In addition, at 2d z = 2 at the unstable 
FP as well. Thus, z^ = z at 2d in general and strong 
dynamic scaling prevails. 


IV. MODEL II: NLS COUPLED WITH A 
STOKESIAN VELOCITY FIELD 

Next, we construct a simple nonequilibrium extension 
of NLS by coupling it to a velocity field v dynamically 
via advection. This may be relevant, e.g., in the dy¬ 
namics of a collection of 2d nematic liquid crystals in a 
2d flow. These are represented by a unit vector in the 
coarse-grained limit [20|, similar to the 0{N) spins (with 
N = 2). Of course, the dynamics of equilibrium or ac¬ 
tive |2l| nematics are far more complex than our simplis¬ 
tic model here; nonetheless, our results in this Section 
may serve as a prototype for the orientational dynam¬ 
ics in systems coupled with an advective velocity field. 
The usual relaxational equation of motion for $ is then 
supplemented by an advective nonlinearity. We obtain 

A T* 

-f Av. V$ = -F^ + e, (47) 

(54) 

where A is a coupling constant. For simplicity we assume 
V to obey a stochastically forced generalized Stokes equa¬ 
tion (neglecting inertia) of the form 

rjV'^Va - Vail = aoVa4> • -f fa. (48) 


in the limit r —?> oo. Thus, the spatial dependence of 
Cs(r) is characterized by a mo del-dependent exponent 
A, which is reminiscent of QLRO in the XY model be¬ 
low the Kosterlitz-Thouless (KT) transition. Nonethe¬ 
less, the spatial dependence in is markedly weaker 
than the algebraic decay of equal-time correlators associ¬ 
ated with QLRO 0. Again for a larger A, g{x) is smaller 
and the spatial decay of C's(r) is accordingly weaker. This 
is consistent with the role of A in favoring order in the 
system. 

We now evaluate the dynamic critical exponent z of 
the broken symmetry modes tt^. The flow equation for 
kF is obtained as 


d dr dn 


( 45 ) 


The first term on the right hand side of (H51) is a 
symmetry-allowed feedback of the 0{N) spin on v. Such 
a feedback may originate from a stress tensor due to the 
0{N) spins Ea /3 = —aoVa^ • such that the stress 

vanishes for $ = const.. This feedback term in (H51) 
is reminiscent of an V-component generalization of the 
feedback term in the generalized Navier Stokes equation 
for a non-critical binary fluid mixture [^ . We allow 
the coupling constant ckq to be both positive or nega¬ 
tive. Further 11 is the generalized pressure. Notice that 
the coupling of 4> with Va is of dynamic origin, unlike in 
Model I, where such couplings are of static or thermo¬ 
dynamic origin (i.e., can be obtained from a free energy 
functional). Model II reduces to pure NLS for v = 0. 
Stochastic force f is assumed to be zero-mean with a 
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Gaussian distribution, having a variance in the Fourier 
space that is given by 

(/a(q,i)//3(-q,0)) = Doq^S{t)5ai3. (49) 

Noise 6 is again chosen to be a Gaussian zero mean white 
noise with a variance given by (0). Note however that 
being out of equilibrium, in (jb)) no longer has the in¬ 
terpretation of the temperature, although D still has the 
same physical dimension as T. Thus, D/k continues to 
be a dimensionless number. In what follows below, we 
generalize the usual low T expansion of equilibrium NTS 
and expand in small D/k. Apart from the formal simi¬ 
larity between D here and T, it is reasonable to expect 
that a low noise, marked by a low value of D, should fa¬ 
vor setting any order; in contrast, a high value of D (i.e., 
high noise) should distablize it. Thus, an expansion in 
D/k is not only a formal generalization of the standard 
low T expansion, it is also physically meaningful. We en¬ 
force incompressibility on v. Thus, 11 may be eliminated 
by imposing V • v = 0 on (|15)) . We eliminate 11 to express 
V as 

= + (50) 

Here, — da9p/^^ is the transverse projection 

operator. Equation (1501) may be used to eliminate v in 
dUl). Evidently, Va enters into the dynamics of $ in two 
distinctly different ways: through (i) the multiplicative 
noise /« and (ii) the deterministic term with the cou¬ 
pling constant ckq. Now proceeding as for Model I and 
replacing v by (l50)) . we obtain the action functional 


correction at the lowest order in D/k. We find 


/DA < 

D D 


(r) 

= ^ + ^A, 

(52) 

/1\< 

1 


(r) 

= p(l + A), 

(53) 

K< 

= /c[l-kA(l-^)], 

(54) 

h< 

= A. 

(55) 


where, A = ^.q^+hi ’ “ Model I and 

/i = Aao/(^??r); superscript < has the same implication 
as in our analysis of Model I above. Since A ^ D/k^ the 
above one-loop corrections are already 0{D/k). Thus, 
we need to find out only 0{D/k)^ corrections, if any, 
to A and ag. We show in Appendix lEl that there are 
indeed no corrections to A and ag at 0{D/ k). To evaluate 
the corrections we once again employ Wilson momentum 
DRG. We scale the fields, g, fl, (pi and pi in the same 
way as for Model I. The parameters are thus scaled in 


the following way: 

N — \ 

/i; = + A], (56) 

(f) = + (57) 

(^) = + (58) 

k' = A(l-/i)]. (59) 

Substituting Eq. (l58l) in Eq. dsn, we obtain, 

D' = iC^Db\ (60) 


Sii = J d’^xdtlY^.^ni - ni{^^ - KV^TTi} + 

1 1 Stt? „ „ 

+ - 2/ii(l-k 7r^/2)}] 

- J d‘‘xdtTri[--^{ypT:j){y‘^Trj)]ya'n-i 

- J d’^xdtXlTi^^fpVaTTj, (51) 


In Model I above, we have argued that 
/m(7fm(q,cj)7rm(-q,-c*;)) and (^(|7rm(q, w)p) must 
scale in the same way due to EDT. We used this to set 
D' = D. Model II does not satisfy EDT. Nonetheless, 
using the freedom to choose one of the rescaling factors 
(^ or ^) freely, we continue to impose D' = D. Thus, 
from Eq. (EQI), ^ Once again arguing that all 

components of hi should scale in the same manner, ^ 
evaluates to b‘^^^[l — A^A]. Hence, ^ = b'^[l — A^A]. 
Using these values of ^ and ^ in Eq. ([59]), we get 


by expanding in (assumed small) D/k, akin to the Model 
I above in Sec. lHIl A = A/T; TTi, i = 1,.., iV—1 is the com¬ 
ponent of the spin fluctuations transverse to the ordering 
direction and is a {N — l)-component vector. Notice that 
there is no analog of of Model I in the present study; 
thus, unlike Model I, we do not restrict ourselves to any 
particular value of D (which is the nonequilibrium analog 
of T here). As before, in order to proceed systematically, 
we use Wilson momentum shell DRG to evaluate the loop 
integrals. Due to the couplings A and ag, there are addi¬ 
tional corrections to k over and above the corrections in 
the pure NLS model. In particular, k receives a O(agA) 


K' = K6'^-2[i_(iV-2)A-^A]. (61) 

We substitute 6 = e*, 1 —>■ 0 and obtain 


^ _ {N-2)D _ fiD 

dl ^ 2TrK 2 ttk 


(62) 


The physical interpretation of (l62ll follows our analysis 
of (l^tl) in Model I closely. As long as N > 2 — p there 
exists a DRG FP for the flow equation (l62l) , such that at 
the FP we have for the reduced noise strength D = D/k 


D* 


27re 

N-2 + n’ 


(63) 
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yielding a (reduced) critical noise strength D* ^ 0{e), 
demarcating between a ferromagnetic phase (ordered 
state) at D < D* and a paramagnetic phase (disordered 
state) at -D > £>* at d = 2 + e. Equation (1551) clearly 

shows (11 = 2 for Model 11. Since the minimum physi¬ 
cally realisable value of is 2, > 2 — p, for ^ > 0 is 

always satisfied and hence a DRG FP always exists. Thus 
for N > 2 — fj,, D* = 0 at 2d, lifting the indeterminacy of 
the transition for A^ = 2 NTS at 2d. What this implies 
about the role of topological defects at 2d for A^ = 2 re¬ 
quires further investigations. Since in general, D* = 0 at 
2d for ^ > 0 and N > 2, there are no ordered phase at 
2d. This is similar to the results for pure NTS at 2d. 

Consider now /r < 0. The DRG FP of (1551) ceases to 
exist for A^ < 2 + \^\. Let us focus at 2d. We define 

a critical by |/ic| = N — 2, < 0, such that at 2d, 

dn/dl = 0 yields an unstable FP given by /i = p-o identi¬ 
cally [see (1551) above]. This again is reminiscent of a sec¬ 
ond order transition at /i = /rd see discussions for Model 
I above. Note that for N = 2, fic = 0. For /r < 
we again obtain a renormalized g-dependent stiffness 
K{q) = — l^jD ln(aog)/(27r) at 2d, as in Model I above. 
Naturally, this implies for the variance Aq = (lnlnL)/|p,| 
in 2d for large L, again similar to Model I. Thus, as /i 
is varied. Model II undergoes a phase transition between 
a soft phase with a vanishing effective stiffness Ke for a 
large enough system and a stiff phase with a diverging He 
in TL. Nonetheless, this transition for Model II is purely 
a nonequilibrium transition, due to the nonequilibrium 
origin of the coupling /r. In contrast, the corresponding 
transition in Model I is an equilibrium phase transition. 
One can analogously define a persistence length C such 
that «((’) = 0, such that C remains finite in the disor¬ 
dered phase but diverges as /x —>• /Xc < 0 at 2d. As in 
Model I, the stiff phase for /x < /Xc at 2d has no true 
LRO; the magnetic order parameter m = {a) vanishes 
for systems in TL both in the ordered and the disordered 
phases at 2d, and hence, cannot be used to describe the 
transition at 2d. Again like Model I, order parameter de¬ 
fined by (1551) may be used to delineate the ordered phase 
from the disordered one. The order parameter exponent, 
unsurprisingly, is 1, same as in Model 1. Notice that for 
p > 0, ( < ({NLS); hence, an increasing p makes the 
system more disordered. The flow equation for the scale- 
dependent reduced noise strength D may be calculated in 
analogy with the flow equation for T in Model I; the as¬ 
sociated physical interpretations for the solutions oi D{1) 
are similar to those for T{1) in Model 1. Specifically at 2d 
for N = 2, 


dD - n iX 

-= D^ — . 

dl 27r 


(64) 


Thus, for any p < 0, D(l) always flows to zero corre¬ 
sponding to the stiff phase, where for any p > 0, D[l) di¬ 
verges as Z —>■ In^, corresponding to the soft phase. This 
is consistent with /Xc = 0 for A = 2. The Debye-Waller 
factor exp(—2IF) and the spin-spin correlation function 


Cs (r) for A = 2 at 2d may be calculated in analogy with 
Model 1. Not surprisingly, W —>• oo at 2d in the stiff 
phase (/X < 0), demonstrating the absence of LRO in the 
stiff phase. In addition, C's(r) ^ [lnr]“^/l^l for ^ < 0 
and r —^ oo at 2d. This is similar to the form of C's(r) 
in Model I stiff phase at 2d, with the mo del-dependent 
exponent being determined by /x < 0. See Fig. [5] for a 
schematic phase diagram of Model II in the N — p plane 
at 2d, highlighting the soft and stiff phases. A corre¬ 
sponding phase diagram of Model II in the e — p plane is 
shown in Fig. |5| 



N 


FIG. 5: (color online)Schematic phase diagram of Model II in 
N — p plane at 2d. The solid inclined (blue) line corresponds 
to the equation N = 2 — p, giving pc. Soft and stiff phases 
are marked (see text). 



FIG. 6: (color online)Schematic phase diagram of Model II in 
e — p plane for N = 2. Symbol, O marks the origin, (0,0). 
Pc = 0 for e = 0. The negative /x-axis (broken) gives the stiff 
phase for e = 0. The inclined line (blue) corresponds to the 
equation p = 2^^ which separates phases with SRO (soft) 
and LRO for /x, e > 0. 
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We now evaluate the dynamical critical exponent, z. 
Using b = e\d = 2 + e and Eq. [551 we obtain the flow 
equation for kF 


d(rK) 

dl 


kT 


z — 2 — e 


{N -2)D 
2ttk 


+ 


kT 


{N - 2)D 
2ttk 


2ttk 


(65) 


At the DRG FP, dK/dl = 0 and hence, therefore, z = 
2 + e — exists. Since for /r > 

D/k = 27re/(7V — 2 + ^), we have z = 2 + 
that depends upon ^ at the DRG FP (similar to the A- 
dependence of z at the DRG FP in Model I). On the 
other hand for /r < /ic, we have z = 2 in the stiff phase. 

In our above analysis, we have treated /i as a constant 
or number. This is justifiable provided it can be shown 
that p. is indeed marginal under the spatial and temporal 
rescaling mentioned above. We now show that below. 
Note that there are no corrections to /i at 0{D/k)^; see 
Appendix|El Thus, ^ is to be affected only trivially under 
the rescaling of space and time. Consider the term in 
Eq. (EJ), A J d’^xdtr^ra'vTTm- In the Fourier space after 
rescaling of momenta, frequencies and fields, this term 
scales as 


- fli - ^ 2 ) X Vi3iq2,^2)iqipTrmiqi,^l)- (66) 

This yields, 

A' = (67) 

where is the corresponding rescaling factor for 
vp{q,V,). Under rescaling q' = bq and U' = we 

can write, 

e^(^(q;,u;)/;3(q',u')) = 2DoqlSiq[U2)mi+^2)b’'^^~^ 

( 68 ) 

yielding where is the rescaling fac¬ 

tor for the noise fa- Now apply the same rescaling of 
wavevectors and frequencies on Eq. (15(71) . yielding in the 
Fourier space 








Padiq') 

m'^ 




q',r2' 

(69) 


We choose ^y = b'^^f = 6 ('^+^+ 2 )/ 2 _ We also define Oq = 
This yields for the scaled coupling 

constant 




/ 


^'Q^O _ f2^-2d-2z 

2r]K'TT 2riKTT 


(70) 


V. MODEL III: NLS COUPLED WITH 
MULTIPLICATIVE NOISES 


Notice that although Model I and Model II above do 
display phase transitions for proper choices of the con¬ 
trol parameters, the magnetic order parameter remains 
zero on both sides of the transition point and there is 
no true LRO in TL at 2d. This sets them apart from 
the usual second order magnetic transitions and raises 
the question what should be a minimal generalized NLS 
with additional degrees of freedom that may show true 
LRO in TL at 2d. True LRO in TL requires finite Aq 
[ see ((551) above] in TL at 2d in the stiff phase. In the stiff 
phases of Model I and Model II, Ke{q) ^ — ln( 5 ao) in the 
long wavelength limit that is responsible for the InlnL 
behavior of Aq. In order to have a finite Aq in TL, Ke{q) 
should diverge more strongly than — hi{qao) in TL at 2d. 
It is reasonable to expect that the presence of long range 
noises in the system may lead to a finite Aq in TL at 2d, 
indicating the existence of true LRO in TL. We study 
this in this Section by using simple reduced models. 

In this Section, we construct simplest possible gener¬ 
alizations of NLS by minimally coupling NLS with mul¬ 
tiplicative noises of given structures. It remains to be 
seen how the explicit forms of the couplings with the 
multiplicative noises affect the long wavelength proper¬ 
ties. Consider first the simple case of a scalar multi¬ 
plicative noise g{x,t), coupled to the 0{N) spin d) via 
symmetry-allowed minimal couplings. The nonconserved 
relaxational dynamics of $ is given by 




(71) 


where A is a coupling constant and 5 is a zero-mean 5- 
correlated Gaussian white noise with a given variance, 
say 1; 0 is a zero-mean Gaussian white noise with a vari¬ 
ance given by ( 0 ). As before we use the parameteriza¬ 
tion $ = (cr, tt) and for an assumed nearly ordered state, 
cr Ri 1. The corresponding generating functional in a 
compact form is 

Zq = j V^V^exp[j d‘*a;dt^$-$ 

+ + (72, 


Now as before, impose d> • d> = 0. This evidently elim¬ 
inates the multiplicative noise term in dZH), leaving it 
identical to that for pure NLS. Thus, the ensuing dy¬ 
namics is identical to that of pure NLS with g having no 
effects on it. 

We now generalize g to vector multiplicative noises and 
study both the nonconserved and conserved dynamics of 
NLS, coupled to vector multiplicative noises Oq, (with d- 
components in a d-dimensional space) via minimal cou¬ 
plings. The corresponding nonconserved and conserved 
equations of motion of NLS are 




establishing the marginality of g. to 0{D/k)^. 


(73) 
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and 

^ + A2V-(a$) =rv2^+V„0„, (74) 

respectively. Stochastic noise 9i is Gaussian-distributed 
with zero mean and variance as given by m, where as 
stochastic noise Qia, again assumed to be zero mean 
Gaussian distributed, has a variance 

(0,„(x,t)0,^(O,O)) = 2Drd{^)Sit)d,,6a,p. (75) 

In addition, Oq, is a vector (of d-components in a el- 
dimensional space) multiplicative noise. We choose Qa 
to be Gaussian-distributed with zero mean and a vari¬ 
ance 

(aa(q,t)a/3(-q,0)) = [mPa/3 +nQa/3]^, (76) 

in the Fourier space. Here, Qap = qaqp/q^ = ^ap — Pap is 
the longitudinal projection operator in the Fourier space, 
Pa /3 is the transverse projection operator defined before. 
Positivity of (|aQ(q,t)p) implies m{d — 1) + n > 0, al¬ 
lowing anyone of m or n to become negative, subject 
to the positivity of (|aa(q,t)P). We choose y, a pos¬ 
itive or negative integer, as a free parameter that de¬ 
fines the scaling of (1751) . The presence of a long-range 
noise coupled to the 0{N) may be realized in a lattice- 
gas type model as follows. Imagine a 2d lattice, where 
sites are occupied by (point) particles, with a spin $ at¬ 
tached to it. Now assume that the particles are noncon- 
served (due, e.g., to evaporation and absorption) - hence 
their density fluctuations are short-lived and drop out 
of the dynamical descriptions of the model in the long- 
wavelength, long-time limit. Now if the particles are sub¬ 
ject to Gaussian-distributed random forces (e.g., if the 
particles are charged and encounter random electromag¬ 
netic pulses, but neglect the inter-particle interactions) 
with appropriately chosen variances, effective dynamical 
equations for the 0{N) spins in this lattice-gas model 
in the long wavelength, long time limit should have the 
form (1731) on symmetry grounds, as the total spin, being 
attached to the particles, is not conserved. On the other 
hand, if there are no particle non-conserving events, but 
the diffusivities of the particles are diverging, then again 
the density fluctuations relax fast (keeping the spins con¬ 
served) and the effective equations of motion for the spin 
variables should be of the structure as in dZl). Notice 
that we do not specify any particular spin-particle in¬ 
teractions above, as our arguments are sufficiently gen¬ 
eral and should hold for any short rang e interactions be¬ 
tween the spins and the particles [^. While this is a 
simplistic motivation and Model III is admittedly arti¬ 
ficially constructed, it is a reduced model that is pur¬ 
posefully designed to study (i) under what conditions, 
true LRO might set in NFS with additional degrees of 
freedom in TL at 2d, and (ii) whether the nonconserved 
and conserved dynamics yield the same result in the long- 
wavelength limit. We will see below that our results from 


Model III are sensitive to the value of y, such that Ag 
may become independent of L leading to LRO in 2d in 
the nonconserved version of Model III, but not in its con¬ 
served version. 

First consider the nonconserved version of Model III. 
The action functional after expanding up to the linear 
order in D/k is (see Si and Sn respectively for Model I 
and Model II above) 

/ JJ 

d‘^xdt—(TTjTTj)^ 

/■ , 1 - r 1 5 2 

+ / d xdt-TTi’nA - 

J 2 " Tdt 

7r2 

-h -2hi{l +—)] 

- Y J d‘^xdtTrm{a-'V)Trm 

P J d’^xdt'ir'^. (77) 

If Ais = Ai/F = 0, one recovers the usual NTS action. 
Non-zero Ais contributes only to D/T to the leading or- 

der at viz., ^^m'^m((a-V)7rn(a.V)7rn). This 

correction evaluates to {n\\g/2) J Hence, 

(f) =|^r("+^)^[l + A + ^A], (78) 

where, A = J ^|^iA/(/cg^). Proceeding identically as in 
Model I and Model H, we obtain for D': 

„\2 p_ 

D' = Db^^C^[l + ^-A]. (79) 

On imposing D' = D, 

„\2 p_ 

e = (80) 

Now, k' = -b A). On substituting ^ = 

A^A], 

t 7A2 r 

k'= nb^--^[l-{N-2)A- -^A]. ( 81 ) 

Notice that there is a contribution to D' with coefficient 
n that originates from the multiplicative noise. Its eval¬ 
uation in details is given in Appendix IfI 

We now have three distinct possibilities, namely y > 
2,y < 2 and y = 2. Clearly, for y > 2, for small q, 
A dominates over A, where as, for y < 2, A is subdomi¬ 
nant. Both terms are equally relevant (in a scaling sense) 
for y = 2. It is important to consider the scaling of the 
parameter y = i that appears as an effective expan¬ 
sion parameter (see below), for an arbitrary y. To find 
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that we look at the scaling of Ais for y ^ 2. Similar to 
Eq. (I5S1) . we can write 

a;, = (82) 


where, is the rescaling factor for the field a. Using 
Eq. (1751) . we evaluate ^a'- 

= yyl%(d+z)i2_ ^ 33 ) 

Using values of ^ and we obtain 

A), = [1 - (iV - 1)A - ^A], (84) 

Thus, up to order (Z1 /k)°, A(g = . Hence, 

A;2r' = XlJ'by-'^ ^ x! = Thus, with b = e' 

and small I 

f = Xta-2). (85) 

Therefore, x is marginal for j/ = 2, but grows (decays) 
indefinitely for y > (<)2. Let us first consider the case 
when y = 2. Using e = d — 2 


dK {N-2 + x)D 

^ = -— 


( 86 ) 


As long asA — 2 + x>0, FP for the flow equation for k 
yields a critical value D* for the reduced noise strength 
D = D/k. We get 


D* 


27re 

A-2 + x’ 


(87) 


demarcating a “high noise” paramagnetic disorder phase 
and a “low noise” ferromagnetically ordered phase at d = 
2 + e. At d = 2, e = 0, D* = 0. This then implies the 
absence of any ferromagnetically ordered state, except 
for at = 0, analog of T = 0 for equilibrium systems. 
This again shows that d^ = 2. 

Notice that negative values are allowed for n, hence 
X can also be negative. Focusing on 2d and with Xc = 
— (A — 2) we write 


large enough system size, allowing us to define a per¬ 
sistence length ( in direct analogy with Model I and 
II above. Similar to Model II, we generally find that 
C(x > 0) < C{NLS) < ((x < 0). The phase transition 
from soft to stiff phase may be described by the same 
order parameter as for Model I and Model II. The flow 
of the reduced noise strength D = D/k may be obtained 
as in Model I and Model II. We find 


dP ^ p2^X 
dl 2 tt ’ 


(90) 


suggesting that for Ay > 0, D flows to zero in TL, im¬ 
plying an ordered state, where as for Ay < 0, oo 

as L -X C, indicating a disordered state. This is sim¬ 
ilar to our analysis for Model I. Notice that Xc = 0 
for A = 2. Then following our analysis for Model I, 
the Debye-Waller factor exp(— 2IU) and the equal-time 
spin-spin correlator C's(r) for A = 2 at 2d may be cal¬ 
culated in the stiff phase at 2d. Not surprisingly and 
similar to the stiff phases in Model I and Model II, we 
find exp(—2IT) —?> 0 for L —oo, precluding any LRO. 
Furthermore, C's(r) ^ l/[lnr]l^l,y < 0 in the stiff phase 
at 2d. Again not surprisingly, the spatial decay of C's(r) 
is the same as in the stiff phases of Model I and Model 
II, and characterized by |y|, a mo del-dependent quan¬ 
tity. A schematic phase diagram of the nonconserved 
version of Model III in the x~P plane is shown in Fig. [7] 
{d = 2,N = 2). 


Stiff 

soft 

X< 0 

X>0 


0 


X 


^ ^ (Xe - x)D ^ Ixcl +X ^ ^ AyP 

dl 2 tt 2tt 2tt ’ 

where Ay = -(lyj -f y) giving 


( 88 ) 


(Ay)D 

K = ---In(qao) -b kq, (89) 

Ztt 

where kq is a microscopic stiffness. For y < 0 and 
|y| > lyd. Ay > 0. Hence, K{q) diverges logarithmi¬ 
cally in TL, g —^ 0. Hence, the variance of the transverse 
spin fluctuations Ag ^ -^InlnL. This is reminiscent 
of our results in the stiff (ordered) phases of Model I 
and Model H, with y playing the role of the control pa¬ 
rameter. For Ay < 0, renormalized k vanishes for a 


FIG. 7: Phase diagram of nonconserved Model III in the y—D 
plane at 2d for A = 2, t/ = 2. Symbol O = (0, 0) is the origin. 
Stiff (y < 0, renormalized D/k 0) and soft with SRO 
(y > 0, renormalized D/k diverges) are marked. 


What happens when y ^ 21 Assume y > 2; the y-term 
in (IFTl) clearly dominates. Noting that j^^^d/^q/qy = 

Z/A^“^, 6 = e* ~ 1 -f /, for y < 0, the recursion relation 
for K now reads 


dK 


= K 


D\xil)\ 


2'!rAy 


(91) 


where e = d — 2. Note that the pure NLS term has been 
dropped in (j9ip since it is less infra-red divergent than 
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the term arriving from the noise for y > 2. Now at 2d 
setting e = 0, dn/dl > 0, implying d^ < 2. Solving 
Eq. (|85ll and Eq. (IMll together we obtain, 


<l) = 


\X\D 


27r(j/ — 2)Ay ^ 






(92) 


where kq is a microscopic stiffness. Therefore, /c(g) di¬ 
verges as a power law in g as g —>■ 0. It should be noted 
that for x > 0 , the recursion relation for k has the form 



Dxjl) 

2TTAy~'^K 


(93) 


Hence the g-dependence of k for % > 0 is given by 


( \ 2-y , 

K(g) = ——q + kq 


(94) 


27r(y - 2) 

Now for y < 0, renormalized variance, Aq is given by 
dq D \x\D dq 1 


An — 


I2tt qK{q) 2 tt 2tt q^ y 

-\ A 


Ixl^ 


27r(j/ - 2) 


1 


72 - 1 / 


(95) 


L-i 


which gives a finite value for Ag in TL, L —>■ 00 , in¬ 
dicating true LRO. Consider the Debye Waller factor 
exp(— 2 IE), as defined above, and the spin-spin corre¬ 
lation C's(r) at 2d for A = 2 with y < 0. Proceeding 
in analogy with Model I, we find that W{= ■^) is finite 
in TL. Thus the order parameter does not vanish in TL, 
indicating LRO in TL. In the same way, 

5 (x) = ^([^(x, t) - 7 f( 0 , t)f) = Ao{L^ 00 ), (96) 


a finite number in the limit r —> 00. Thus, Cs(r —>■ 00) 
is finite, consistent with a finite W as obtained above. 
This firmly establishes LRO (ferromagnetic order) for 
X < 0 , 2 / > 2 at 2d. Notice that for any y > 2 and 
X < 0, the system is always in the ordered phase; there 
is no phase transition at any d. This is because renor¬ 
malized D/k —!► 0 in TL at any d for 2 / > 2, x < 0. The 
existence of LRO should not come as a surprise, since 
nonconserved Model III is driven out of equilibrium by 
the long-range multiplicative noise. We now briefly con¬ 
sider the consequence of x > 0. Clearly, from (IMl) . k{1) 
rapidly drops to zero as the system size approaches the 
persistence length given by 


c 


Kpjy - 2 ) 
xD 


1 / 0 - 2 ) 


(97) 


Thus, nonconserved Model III with y > 2 forms an ex¬ 
ample of dynamic second order phase transition between 
a disordered state with SRO and an ordered state with 


LRO with X as a tuning parameter. In dynamic phase 
transitions, a model goes from one phase to another 
upon tuning a model parameter of dynamic origin. Well- 
known examples of dynamic phase transitions include the 
Kardar-Parisi-Zhang (KPZ) equation for surface growth 
at dimensions d > 2 [ 2 ^ and active-to-absorbing state 
phase transitions . Our nonconserved Model III with 
2 / > 2 is an analog of these examples where a continuous 
symmetry is dynamically broken. 

Lastly, for 2 / < 2, the pure NLS contribution in (ISTl) 
dominates in the long wavelength limit, regardless of the 
sign of X- Thus, in that limit the properties of noncon¬ 
served Model III with 2 / < 2 is identical to that of pure 
NLS, with no order at a finite T in 2d. A phase diagram 
of the nonconserved version of Model III in the x~y plane 
is shown in Fig. [5] with d = 2, N = 2. Below in Fig. 
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FIG. 8 : Phase diagram of nonconserved Model III in x — 2 / 
plane at 2d with N = 2. Symbol O denotes the origin (0,0). 
LRO, soft and stiff phases are marked. 


we provide a phase diagram of nonconserved Model III 
in the y — e plane with x < 0 - 

We now calculate the dynamic exponent z at the DRG 
FP and in the stiff/LRO phases. For y < 2, results of 
the pure NLS ensue, precluding any stiff phase at 2d for 
all A > 2. On the other hand, for 2 / > 2, x < 0 we have 
for the recursion relation for P 


P'= _ |;^|A]. (98) 


The recursion relation for P for 2 / > 2, thus takes the 
following form: 


\x\D{i-by-^) - 

2tt{ 2 — y)Ay~‘^K 


(99) 


As before, the dynamic exponent z is obtained from the 
DRG flow equation for kP. Considering y > 2, using (mD 
and proceeding as for Model I and Model II, we find z = 2 
in the LRO phase. Precisely at y = 2, again proceeding 
as for Model I and Model II, we have for x < 0 


z — 2 -f e — 


(A-2 -|x|)D 

2ttk 


( 100 ) 
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The contribution to from the multiplicative noise 
may still be represented diagrammatically as in Ap¬ 
pendix If] However, its explicit expression is clearly dif¬ 
ferent from the corresponding form for the nonconserved 
Model III. Using the same scaling procedure for q, tt 
and TT as above, we can now write ^ 

These can be used to write the following re-scaled cor¬ 
rected values of the parameters: 





- A- 


Ai 


2r£i d‘ 


(106) 




FIG. 9: (color online)Phase diagram of nonconserved Model 
III iny — e plane. Regions in the phase space with LRO, SRO 
(soft phase) and stiff phase (circle) are marked. 


at the unstable DRG FP for |x| < \xc\- However, for 
Ixl > \Xc\, z = 2 in the stiff phase, as in Model I and 
Model H. 

We now study the conserved version of Model HI. The 
corresponding action, Snicon is given by 

Silicon = J 

- J d’^xdtitil-Y^^dtn^ - KV’^TTi] 


+ J d’^xdt^iViTTjTTj))^ 

f ,d A r 1 d 2 

+ / d xdt-^,.A^-. 


n' = -h A]. 


(108) 


Using Eq. (11061) and Eq. (11071) and imposing D' = D as 
before, we obtain 


^=b^i[i+p,^], 


(109) 


where. Pi = Now using the values of ^ and 


k' = k5'^-2[i_(7V_2)A-Pi—]. (110) 

Using values of ^ and we obtain 
(h). = hi,^=4“li_,jv-i)A-p4| ,111) 


- 


A 2 

T 


J d^xdt^^{a ■'V)Ttm + p J d’^xdtn'^ (lOlf 


2h (1 -I- up to order {D/k)°, (^)'^ = ^+y+‘^_ Hence, 

2 ^ (M)2r5y p/ = p^i,y .,^ith 6 = e' 


r ; ^ t r 
nd small I 


The leading contributions to D/T come from (i) 
ai'^jT^j) and (ii) a correction second or¬ 
der in Aa/F, viz., 
pendix O 

For the conserved model, the correlator, (niiTj) = 

r d'^q dO. 2DSij _ r d’^q Dr.. _ A r. . p 

J (27r)<i 271- r 52 (_u^+„ 2 '^ J (2ii-p D V 

small hi). Furthermore, {da'^iap'Kj) evaluates to 


dPi 

~dr 


= yPi 


Pi{l) = Piexp[yl]. 


( 112 ) 


I 


d^q D 


■[m{d — 1) -I- n] 




= A- 


values, we find 

(f 


-. Using these 


Therefore, Pi is marginal for y = 0, but grows (decays) 
indefinitely for y > (<)2. 

The second and third terms in the RHS of (jllOl) com¬ 
pete when y = Q. For ?/ < 0, the pure NTS behavior 
follows while for y > 0, the nonequilibrium contribution 
dominates at small q. For d = 2, 

where, b = = 1 + 1 and A is an upper wavenumber 

cut-off. Similar to the non-conserved case (see above), 
here we find for y > 0 


r 

A Ai A, 

= 1-bA-f ^ , 

F ^ 2FP d^’ 

(102) 

1 

f< 

= ^[1 + A], 

(103) 

«< 

= k[1 + A], 

(104) 

h< 

N -1 

= hi[l + ^-A]. 

(105) 


dn 

U 


= K 


e-P 


D 


27rA*^/c 


(113) 


where e = d — 2 and P = lniilLJl+nl p^ _ Hence dP/dl = 
[mid-l)+n] dp^/dl = yp At d = 2, 


dn 

m 


PD 


2t: ky' 


( 114 ) 
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Solving the flow equations for P and k together, 


VI. SUMMARY AND OUTLOOK 


<l) = - 


PD 

2TTyhy 


PD 


—9 ^ 

2TTy 


Kq. 


(115) 


Thus as 9 —>■ 0, renormalised K{q) runs away rapidly to 
large negative values, implying the absence of any stiff 
phase in TL. The persistence length, C at which K{q) 
vanishes is given by 


c 


'Koy 


i/y 


IPDi 


(116) 


decreasing rapidly as y increases for fixed ko,P and D\ 
thus n{l) rapidly vanishes to 0 at scale ^ C,. Clearly, 
C < C{NLS) for a large enough kq. For y < 0, evidently 
the results of the pure NTS holds. A schematic phase 
diagram of the conserved version of Model III in the y—P 
plane is shown in Fig. fTOl below. 
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FIG. 10: Phase diagram of nonconserved Model III in y — P 
plane at 2d . Symbol, O = (0, 0) denotes the origin. For both 
positive and negative values of y, the system remains soft but 
with different values of as marked in the figure. 

The take home message from this Section is that 
merely the presence of long-range noise does not auto¬ 
matically ensure phase transition to LRO in the system. 
The detailed mechanism of coupling of the noise with 
the 0{N) and the underlying dynamics (conserved versus 
nonconserved) as well as the structure of the noise (e.g., 
scalar versus vector) are also equally important. Signif¬ 
icant differences between non-conserved and conserved 
dynamics of Model III are not unexpected, as for sys¬ 
tems out of equilibrium, the static or equal-time proper¬ 
ties obtained from conserved and nonconserved dynamics 
are not necessarily same, i.e., the equal-time properties 
depend on the underlying dynamics. This is in contrast 
to the well-known equivalence of ensembles in equilibrium 
systems. 


In this article, we have addressed phase transitions and 
the nature of order, if any, in certain generalized NTS 
constructed by us, where the 0{N) spins are minimally 
coupled to additional degrees of freedom. To this end, 
we set up the dynamics of the nonconserved NTS with 
(i) thermodynamically coupled Ising spins (Model I), (ii) 
a dynamically coupled Stokesian velocity field (Model 
II) and (hi) a dynamically coupled vector multiplicative 
noise (Model III). We also study the conserved version 
of Model III. In Model I, the couplings are of thermody¬ 
namic origin, where as in Models II and III, the couplings 
with the additional degrees are of dynamical origin. We 
use a low noise variance expansion^ that for Model I is 
identical to the well-known low-T expansion for the pure 
NTS, where as for Model II and Model III, it is a gener¬ 
alization of the standard low-T expansion. Our studies 
uncover unusual phase transitions in the models as the 
relevant couplings exceed particular threshold or critical 
values; in Models I and II, and in the nonconserved ver¬ 
sion of Model III for ?/ = 2, the transition is between 
a soft phase with SRO, in which the renormalized spin 
stiffness n vanishes in the long wavelength limit and a 
stiff, ordered phase, where k diverges as InT in 2d. In 
the latter phase, the variance Aq of the transverse spin 
fluctuations scale proportional to In InT with a model- 
dependent proportionality constant at 2d, a dependence 
of L that is substantially weaker than the InT depen¬ 
dence of the broken symmetry modes in classical elastic 
Hamiltonians that exhibit QLRO. This implies a spatial 
decay of the equal-time spin-spin correlation function in 
powers of inverse logarithm of the (large) spatial sepa¬ 
ration, characterized by model-dependent exponents, at 
2d for N = 2. This mo del-dependence of the exponents 
are reminiscent of the mo del-dependent exponents that 
characterize the algebraic decay of the spin-spin equal¬ 
time correlator in the equilibrium XY model displaying 
QLRO; nevertheless, these power laws of inverse loga¬ 
rithm are distinctly slower varying function of distance 
than the algebraic decay in QLRO. Notice that this tran¬ 
sition is not temperature driven, but rather controlled by 
the coupling constants at fixed T or at the fixed noise am¬ 
plitudes. These transitions have no analog in pure NLS. 
Nonetheless, the magnetic order parameter remains zero 
even in the stiff phase in TL at 2d. For y > 2, Aq becomes 
independent of T even at 2d. This implies LRO in the sys¬ 
tem. The presence of the long-range multiplicative noise 
makes the model nonequilibrium and takes it out of the 
validity of MWT. With y as the tuning parameter, our 
nonconserved Model III with y > 2 displays a dynamic 
second order phase transition. We emphasize that there 
is indeed no clash between our results and MWT. Notice 
however, the mere presence of long-range noise itself does 
not necessarily lead to LRO, e.g., a scalar multiplicative 
noise minimally coupled to the 0{N) spins have no effect 
and the model is identical to pure NLS. In addition, the 
conserved version of Model III does not display any LRO 
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for y > 0. The lesson, therefore, is that the specific mi¬ 
croscopic couplings between the noise and the 0{N) spins 
may lead to LRO. Furthermore, the conserved version of 
Model III does not display any phase transition at 2d, 
unlike its nonconserved counterpart, and has long wave¬ 
length properties that are qualitatively similar to pure 
NTS. This highlights the lack of equivalence between the 
conserved and nonconserved dynamics of Model III, em¬ 
phasizing the underlying nonequilibrium nature of the 
dynamics. We also obtained the dynamic exponent z at 
the DRG FP and in the stiff phases of Model I, Model 
II and nonconserved Model III, and in the LRO phase of 
nonconserved Model III with long-range noise. 

Our models here are admittedly artificially constructed 
and simple. Nevertheless, they reveal important physical 
insight for more realistic but complicated models. We 
note that in general, the additional degrees of freedom 
can affect the dynamics of the 0{N) spins in two dif¬ 
ferent ways: (i) through new effective deterministic cou¬ 
plings for the 0{N) spins, and (ii) modification of the 
noises, e.g., generation of effective noises in the 0{N) 
spin dynamics. These features should be generic in NTS 
coupled to other fields. Thus, the basic features of our 
results should hold. While the Ising spins in Model I 
have independent dynamics, the overall structure of the 
coupled dynamical equations is constrained by the FDT. 
Instead in Model II, the velocity field being assumed to 
be Stokesian, has no independent dynamics, or, in Model 
III, the multiplicative noise is entirely prescribed by its 
variance. An obvious generalization would be to con¬ 
sider nonequilibrium coupling with additional fields hav¬ 
ing their own dynamical evolutions, e.g., coupling with 
a growing surface described by a generalized KPZ equa¬ 
tion [^ . In this model, whether or not one would find 
true LRO or an analog of the stiff phase, cannot be im¬ 
mediately predicted. Lastly, we have entirely ignored the 
topological defects. It is known that in 2d, such defects 
lead to the Kosterlitz-Thouless transition in the equilib¬ 
rium XY model, where the system undergoes a transition 
from low-T QLRO with bound pairs of defects to high-T 
SRO with free defects @ . It would be interesting to find 
how the results of the present study may change if the de¬ 
fects are accounted for in the models used here. We look 
forward to detailed theoretical studies in this direction. 
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Appendix A: Action for Model I 

We begin with the action functional Sj for Model 1. 
We perform the following rescalings: (i) rescale t —>■ t/n, 
(ii) absorb a factor D/k in ifm, (in) absorb a factor of 
y/n/D in tt^, (iv) tp —>■ '/D'tp,{v)'tp —>■ ■ipK/'/I). This en¬ 
sures that the equal-time correlators (7rm(q, t)7rm(—q, t)) 
and are independent of the model pa¬ 

rameters, and hence 0{D/k)'^. We find for the rescaled 
action 

/ 2 I 

d xdt^ — TTm'ii'm ffml —cliTTm V TTm 

“2X0^ ^ , 2 ^ ID,, ,2 

tv L Ki 

I 1 » / 1 o 2 I v72 2 ^1 2\ 

2 K I K 

hi , /" , 1 -D 2 

— - TTruTTm + d xdt -pTT 

K J K K 

-f ( d‘^xdt[ipip^+ rip 

J I 2 ^2 

where F 2 = ^2 = 1/k. Clearly, all the nonlinear 

terms are 0{D/k). This justifies our perturbative ex¬ 
pansion in small D/k. Thus, the bare propagator and 
correlators, which contain only terms that are 0{T/k)^, 
after Fourier transforming in space and time are 


{tT^TTj) = 

5ij 

(A2) 

=^ + Kq^ + hi' 

II 

2D/r y 

(A3) 

^ + {Kq^ + hi)^ 

(iplp) = 

1 

(A4) 

—iui + r + q^^ 

{ 1 P 1 P) = 

2T 

(A5) 

uj^ + [r + q^r 


Appendix B: Action for Model II 

Use the same scaling for and TTm as in Sec. lAl above 
to obtain 

d xdt\^ —77.^171 jYi ^ '^m 

2XD^,,.,^ ID,, ^2 

V ,3 \1p V f T y., j 

tv L tv 

, , I r, 2 ^2 2 

+ ■7:—T^mT^ni{-7-dtTT -f kV TT 

2 tv i 

- (ft.i/K)7r^) - (ft.i/K)7r™7r^] 

-I- [ d^xdt - pTT^ 

J K K 

- Xao J d'^xdt{D/K)[TTi^^{yi3Trj){V'^Trj)yaT^i], 

(Bl) 
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Appendix C: Action for Model III 

Consider the nonconserved Model III with a vector 
multiplicative noise a. Averaging the generating func¬ 
tional over a leads to a term in Sjunon of the form 

Ajg J dfxdtdf'x',t') 

X ((a(x,t) • V)7rm(a(x',t') • V)7r„(x',t'))a, (Cl) 


respectively. In each of the above cases, renormalization 
of K may be calculated to the leading order in H/k, which 
are identical to the results obtained for Model I, Model 
II and Model III above with N = 2. Equal-time spin 
correlation functions C's(r) as defined above may be cal¬ 
culated in straight forward ways yielding same results as 
above. 


where (...)a implies averaging over the distribution of a. 
Noting that a is J-correlated in time and then using the 
same scaling as in Sec.lBlabove. the new term (EU above 
gets a scale factor {\\g/D){D/ k). Averaging over a yields 
additional factors containing m and n; see (US]). This 
yields a factor ^ for the contribution to D. All 

other terms in Smnon are same as those in the action 
for pure NTS. This establishes the expansion of Smnon 
to 0{D/k). The action Smcon for conserved Model III 
may similarly be expanded up to 0{D/k). 

The correlator and propagator terms for conserved 
Model III are given by: 


< TTniT^n ^ — 


2D/T 




+ + hiY 

dmn 


— l<jJ 


-I- -I- hi 


(C2) 

(C3) 


Appendix D: Alternative formulation for N = 2 

In case of 0(2) spins, using = 1, we write $ = 
exp[i4']. Here, phase d' is a Goldstone or broken sym¬ 
metry mode. Any variation in $ is incorporated as a 
variation in the phase d/. The equations of motion in 
terms of 'k can thus be written as : 

Model I : 


Appendix E: Renormalization of k in Model II 

The term f d'^xdtAn^^f^Va'^'i, in the second order 
evaluates to zero and hence does not contribute to -p-. 

We now show that there is no renormalization to A at 
the one-loop order in Fig. [TTJ 



FIG. 11: Feynman diagram for one loop correction to A in 
Model II. 

This evaluates to 


— =rlKV^'I>-2AV(^^Vi')] + e^. (Dl) 

Model II : 

-f Av.V^-= TkV^T-(-0^,. (D2) 

Model III : Non conserved case 

Aia.V^' = rKV^T-f 0VI,. (D3) 

Model III : Conserved case 


d'^qdVLPapiqjPapiq) 


(/7(q: ^^)/5(-q,-w)) 


2 4 
rj^q^ 


X gi^(7r;,(-q, -n)Tr^{q, H))(7r,.(-q, -H)7r^(q, H)) 

(El) 


which vanishes due to causality. Now, we show that the 
coupling constant ao has no correction to 0{D/k)^. The 
relevant diagram is shown in Fig. dUD. 

Fig. dEI) gives the following contribution to k 


- 


d‘^qidnTr,Pap{q_i){\TTj{q -qi,uj- H)p), (E2) 


— + A2V.(a^') = rv2[-K(v2^')] -p V„0^„///. (D4) 

Here, stochastic noises and 0,j,i are zero-mean Gaus¬ 
sian white noises related to the noises 6 and 0^ above 


as available in Section ITVl above. 

Clearly, the above diagram being proportional to 
{tToiT^p) is proportional to D/k. Thus, there are no cor¬ 
rections to 0{D/k)^. 
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FIG. 13: Feynman diagram for the relevant correction in the 
non-conserved version of Model III. 


The one-loop diagram for the contributions upto 
0{D/k) to -p in the non-conserved version of Model III 


Appendix F: Contribution of the multiplicative 
noises to one-loop corrections for D in Model III 
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